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We calculate analytically and numerically the Dyer-Roeder distance in perfect fluid quintessence 
models and give an accurate fit to the numerical solutions for all the values of the density parameter 
and the quintessence equation of state. Then we apply our solutions to the estimation of Hq from 
. multiple image time delays and find that the inclusion of quintessence modifies sensibly the likelihood 

distribution of Hq , generally reducing the best estimate with respect to a pure cosmological constant. 
£S) , Marginalizing over the other parameters (O to and the quintessence equation of state) , we obtain 

Hq = 71 ± 6 km/sec/Mpc for an empty beam and Ho = 64 ± 4 km/sec/Mpc for a filled beam. We 
also discuss the future prospects for distinguishing quintessence from a cosmological constant with 
time delays. 
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I. INTRODUCTION 
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■ Quintessence (Caldwell et al. 1998) or dark energy is a new component of the cosmic medium that has been 
introduced in order to explain the dimming of distant SNIa (Riess et al. 1998; Perlmutter et al. 1999) through an 
t-H ■ accelerated expansion while at the same time saving the inflationary prediction of a flat universe. The recent measures 
of the CMB at high resolution (Lange et al. 2000, de Bernardis et al. 2000, Balbi et al. 2000) have added to the 
motivations for a conspicuous fraction of unclustered dark energy with negative pressure. In its simplest formulation 
(see e.g. Silveira & Waga 1997), the quintessence component can be modeled as a perfect fluid with equation of state 
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with to in the range 0<to<3(to<2 for acceleration). When m = we have pure cosmological constant, while for 
+-> ■ to = 3 we reduce to the ordinary pressureless matter. The case to = 2 mimicks a universe filled with cosmic strings 
(see e.g. Vilenkin 1984). More realistic models possess an effective equation of state that changes with time, and can 
be modeled by scalar fields (Ratra & Peebles 1988, Wetterich 1995, Frieman et al. 1995, Ferreira and Joyce 1998), 
possibly with coupling to gravity or matter (Baccigalupi, Perrotta & Matarrese 2000, Amendola 2000). 
, The introduction of the new component modifies the universe expansion and introduces at least a new parameter, 
to, in cosmology. Most deep cosmological tests, from large scale structure to CMB, from lensing to deep counting, 
are affected in some way by the presence of the new field. Here we study how a perfect fluid quintessence affects the 
Dyer-Roeder (DR) distance, a necessary tool for all lensing studies (Dyer & Roeder 1972, 1974). The assumption of 
constant equation of state is at least partially justified by the relatively narrow range of redshift we are considering, 
z G (0.4,2.5). We rederive the DR equation in quintessence cosmology, we solve it analytically whenever possible, 
and give a very accurate analytical fit to its numerical solution. Finally, we apply the DR solutions to a likelihood 
determination of Ho through the observations of time-delays in multiple images. The dataset we use is composed 
of only six time-delays, and does not allow to test directly for quintessence; however, we will show that inclusion of 
such cosmologies may have an important impact on the determination of Hq with this method. For instance, we find 
that Hq is smaller that for a pure cosmological constant. In this work we confine ourselves to flat space and extremal 
values of the beam parameter a; in a paper in preparation we extend to curved spaces and general a. 



II. THE DYER-ROEDER DISTANCE IN QUINTESSENCE COSMOLOGY 

In this section we derive the DR distance in quintessence cosmology, find its analytical solutions, when possible, 
and its asymptotic solutions. Finally, we give a very accurate analytical fit to the general numerical solutions as a 
function of Q m and to. 

First of all, let us notice that when quintessence is present, the Friedmann equation (the 0,0 component of the 
Einstein equations in a flat FRW metric) becomes (in units G — c = 1) 
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where Hq is the present value of the Hubble constant, U, m the present value of the matter density parameter, and 
where the scale factor is normalized to unity today. In terms of the redshift z we can write 



where 
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The Ricci focalization equation in a conformally flat metric (such as the FRW metric) with curvature tensor R a p 
is (see e.g. Schneider, Falco & Ehlers 1992) 



y/A+^R af }k a k VA = 
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where A is the beam area and k a = is the tangent vector to the surface of propagation of the light ray, 
and the dot means derivation with respect to the affine parameter v. Multiplying the Einstein's gravitational field 
equation i? Q ^ — ^g a pR = %irT a p by k a k^ and imposing the condition g a pk a k^ = for the null geodesic we obtain 
R aj3 k a k^ = SirT a pk a k^ ; from @j we obtain 



VI +ATiT af} k a k f) \ r A = 0. 
Considering only ordinary pressureless matter and quintessence the energy-momentum tensor writes 

T a fi = p m u a U{) + p q (1 + 7) u a up - rypqgap; 
multiplying by k a k@ , putting 7 = y — 1 and inserting in (|^) we have 
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Now, the angular diameter distance D is defined as the ratio between the diameter of an object and its angular 
diameter. We have then D oc \J~A. Since u a k a = u> = -4- (1 + z), and defining the dimensionless distance r = DHq , 
Eq. (0) writes 



-Pq 



r = 



where we introduced the affine parameter a defined implicitcly by the relation 

g(z) = — = (l + z) 2 E(z) 
where E (z) is defined in Eq. (||). Finally we get the DR equation with the redshift as independent variable 
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The constant a in Eq. ( JlOj ) is the fraction of matter homogeneously distributed inside the beam: when a — all 
the matter is clustered (empty beam), while for a = 1 the matter is spread homogeneously and we recover the usual 
angular diameter distance (filled beam). Notice that in our case the "empty beam" is actually filled uniformly with 
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quintessence. Since the actual value of a is unknown (see however Barber et al. 2000, who argue in favor of a near 
unity), we will adopt in the following the two extremal values a — and a = 1. The appropriate boundary conditions 
are (see e.g. Schneider, Falco & Ehlers 1992) 
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Defining r (0, z) = r (z) these become 

' r (0) = 

Analytical solutions. Eq. ( fTo| ) has an analytical solution only for some values of a and m. Here we list all the 
known analytical solutions. All the cases for a — and m ^ or Sl m ^ 0, 1 are new solutions. For the case m = 
see also Demiansky et al. (2000). 

Case a = 0, m = 

When 77i — the quintessence fluid is the cosmological constant; the DR equation becomes 
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whose solution in integral form is 
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where F [a, 6; c, x] is the Gauss hypergeometric function. 

Case a — 0, m = 1 
The solution is 
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Case a = 0, m = 2 

The solution can be written in terms of the Meijer G function, but the expression is so complicated that it is not 
worth reporting it here. 
Case a = 0, m = 3 

Now quintessence coincides with ordinary matter; however, the choice a = implies that the ordinary matter is 
completely clustered, while quintessence remains homogeneous. The solution is 
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For fi TO = the DR equation is 



(l + ^) : 



d 2 r 777 + 4 



dz 2 



, dr m 
2 (1 + z) ^ + y = ° 



(19) 



3 



and its solution is 



r (z) = < 
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while for fl m = 1 we get 
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The case a = 1 is actually the standard angular diameter distance in a homogeneous universe. We list some particular 
solution here for completeness (see also Bloomficld- Torres & Waga 1996, Waga & Miceli 1999). 

Case a = l,m^3 
The general solution is 



r(z) 



which gives for any m 
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Case a = 1, m = 3 
The solution reduces to 



Case a = 1, TO = 0; 1 
For f2 m = we get 



r (z) = 2 
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(identical to the case a = 0, f2 m = 0). For TO = 1 we reduce to the case a = 1, m = 3. 

Asymptotic limits. The DR equation can be solved analytically in the limit of small or large z. Here we give 
these limits; they will be used to produce accurate analytical fits to the numerical solutions. For small z and a = 
we have 
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whose solution is 
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or, to second order in z 
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For completeness, we quote also the result in the case a ^ 
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The limit to second order is identical to Eq. (p0|) , which shows that for small redshift the degree of emptiness is not 
relevant. 

For large z (and m ^ 3, f2 m ^ 0, 1) the DR equation for a — reduces to 
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If we define r; (z) = ^^^5/4 we obtain the equation 
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whose solution can be written in terms of Bessel functions of first kind: 
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Notice that the large- z limit is always a constant: 

lim n{z) = C 2 

z — >oo 

Since in the limit of large z the C\ term in (^5|) is negligible, we will consider in the following only the Ci term. 
Again for completeness we quote the analogous result for any a^O: 
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Numerical fits. The rest of this section focuses on the a = case, for which we do not have a closed solution. 
Although in the application of the next section we employ the exact numerical solutions, it may prove practical to 
produce analytical fits. We use the asymptotic limits above to build an analytical fit to the full numerical solution, 
valid for all f2 m and all m. We search a fit in the form 



r ftt (z) = r s (z) T s (z) + n (z) Ti (z) 



(37) 
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where the function T s (z) is a step- like curve chosen so as to interpolate from 1 to as z goes from to infinity, and 
Ti(z) interpolates similarly from to 1. We chose 



T s {z) = \ 
Ti (z) = \ 



l-tanh^ln^)) 
l + tanh(iln(^)) 



(38) 



so that the transition occurs at zq and A sets its steepness. We have then three parameters (C2,z ,A) to fit as 
functions of Q m and m. The result for Ci in terms of a simple polynomial fit is 



C 2 = -0.034rj^m 2 - 0.304^m + 0.112O m m 2 + 0.339r4 
-0.078m 2 + 0.397ft m m - 0.6ft m - 0.117m + 0.682. 

For z and A we find convenient to use the following functional form: 

(a$m 5 + aim 4 + a 3 m 3 + a 2 m 2 + aim + a ) sinh (^O 2 ^ + &if2 m 

The values of the fit parameters for z and A are listed in Table I. 

Table I 



(39) 



(40) 





a 5 


a 4 


a 3 


«2 


a x 


a 


b 2 




bo 


zo 


-0.340 


1.775 


-3.324 


2.692 


-0.884 


0.132 


2.005 


-3.025 


4.935 


A 


-0.280 


1.687 


-3.747 


3.725 


-1.564 


0.367 


-0.521 


0.756 


1.998 



Such fits are accurate to better than 5% over the range m e (0, 2) and Q m e (0, 1) as can be seen from Fig. 2, where 
we compare our fits to a sample of exact numerical solutions. 

As a cautionary remark, let us notice that the assumption of a constant a over a very large range in redshift is 
certainly problematic. The results of the next section, however, are obtained in a relatively narrow range of redshifts, 
so that the approximation should be acceptable. 



III. MEASURING Ho IN QUINTESSENCE COSMOLOGY THROUGH TIME-DELAYS 

Now that we are in possess of the general angular diameter distance in quintessence cosmology we can apply it to 
real data. As a first application, in this section we use six observed time-delays to measure H taking into account 
the presence of quintessence fields. 

Let us first present the data. There are only seven lens systems with measured time-delays: B0218+357 (Biggs 
et al. 1999; Lehar et al. 1999; Patnaik Porcas & Browne 1995), Q0957+561 (Bar-Kana et al. 1999; Kundic et al. 
1997), HE1104-1805 (Courbin, Lidman & Magain 1998; Wisotzki, Wucknitz, Lopez & S0rensen 1998), PG1115+080 
(Bar-Kana 1997; Impey et al. 1998; Schechter et al. 1996), B1600+434 (Burud et al. 2000; Koopmans, de Bruyn, 
Xanthopoulos & Fassnacht 2000), B1608+656 (Fassnacht C. D. et al. 1999; Koopmans & Fassnacht 1999) and 
PKS1830-211 (Lehar et al. 1999; Lovell et al. 1998; Wiklind & Combes 1999). Due to the image multiplicity, we have 
in total ten time-delays. The lens model we use, the isothermal lens of Mao, Witt and Keeton (Mao, Witt & Keeton 
2000) cannot be adapted to Q0957+561 and B1608+656 so we are left with five lens systems and six time-delays, as 
detailed in Table II. 



Table II 



Lens / components 


Zd 


Zs 




* (") 




f ( Mpc \ 
J (km/s J 


A/ 
/ 


B0218+357/BA 


0.68 


0.96 


(10.5±0.2)d 


0.24±0.06 


0.10±0.06 


0.031 


87% 


HE1104-1805/AB 


0.77 


2.32 


(0.7±0.1)yr 


2.095±0.008 


1.099±0.004 


0.0108 


16% 


PG1115+080/AB 


0.31 


1.72 


(11.7±1.2)d 


1.147±0.025 


0.950±0.004 


0.0051 


25% 


PG1115+080/CB 


0.31 


1.72 


(25.0±1.6)d 


1.397±0.004 


0.950±0.004 


0.00433 


8% 


B1600+434/BA 


0.42 


1.59 


(47±6)d 


1.14±0.05 


0.25±0.05 


0.0064 


23% 


PKS1830-211/BA 


0.89 


2.51 


(26±5)d 


0.67±0.08 


0.32±0.08 


0.0095 


64% 
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Adopting the isothermal model, the relation between the distances of the deflector (subscript d) and of the source 
(s) and the time-delay between two images labelled 1 and 2 in terms of observable quantities is 

M= l + ZdDdDs M\ (41) 
2 Dds 

where 

At? 2 = [d\ - (42) 
i?i ; 2 is the angular distance between one of the two images and the center of the deflector and where 

D a ,b = -^-r(z a ,z b ;n m ,m). (43) 
Ho 

To estimate H we build the likelihood function 

(f iit (H ,Q m ,m) - f^of 



21og L(H ,n m ,m) = ^2 



2af 



(44) 



where we separated between quantities that contain theoretical parameters and purely observational quantities by 
definining the variables 

k, - • (45) 



ds 



'--{ik^l- (46) 

and where <7j is the total error on /j ;D , obtained by standard error propagation (this error is dominated by the 
uncertainty on the angular positions $1,2)- Then we marginalize over Cl m ,m (defined in the range 0,1 and 0,3, 
respectively) and produce the marginalized and normalized likelihood 

L(H ) = N [ dSl m f dmL(H ,n m ,m), (47) 
Jo Jo 

that represents the likelihood of H (normalized to unity by the constant N) given any possible perfect fluid 
quintessence model. We also evaluated L(m) and L(fl m ), marginalizing over the other variables, but the likeli- 
hoods we obtain are too flat to derive any interesting conclusion on the quintessence parameters. We also estimated 
the effect of imposing a Gaussian prior on f2 m , with mean 0.3 and standard deviation 0.1. We compared the marginal- 
ized likelihood with the likelihood in the "standard" cases Cl m = 0.3 and m = (pure cosmological constant), and 
TO = 1 and m = 3 (no quintessence). Table III summarizes the likelihood results: here l,2,3er stand for a probability 
of 68,95,99%, respectively; the first four cases are marginalized over O m and to, and the prior is the above mentioned 
Gaussian on Cl m . 

The strongest effects are obtained in the empty beam case, a = 0, because in this case the quintessence term is 
dominating in the last term of the DR equation. The likelihood is shifted to lower values, in comparison to the two 
"standard" cases (see Fig. 5), with an increase in the variance. For a = 1 the shift is less evident and there is a 
degeneracy between marginalized likelihood and no quintessence likelihood. Thus the dependency of H for TO and 
to is strongest in the empty beam case. With no prior we obtain 

H = 71 ± 6 km/sec/Mpc if a = 

H a = 64 ± 4 km/sec/Mpc if a = 1 (48) 



Table III 



Model 


Max 


la 


2a 


3a 


a = 0, margin., no prior 


71 


65<H < 77 


60<H < 83 


57 < H < 88 


a = 0, margin., prior 


69 


64 < H < 74 


60 < H < 80 


57 < H < 85 


a = 1, mar gin., no prior 


64 


60 < H < 69 


56 < H < 74 


54 < H < 78 


a = 1, margin., prior 


64 


60 < Ho < 69 


57<H < 74 


54 < H < 78 


a = 0, n rn = 0.3, to = 


73 


69 < H < 78 


65 < H < 84 


62 < H < 89 


a = 0, Q m = 1, to = 3 


75 


70<H < 80 


66 < H < 86 


64 < H < 91 


a = 1, il m = 0.3, to = 


67 


62<H < 72 


59 < H < 77 


57 < H < 81 


a = 1, fl m — 1, to = 3 


63 


59<H < 68 


56<H < 73 


54 < Ho < 77 
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Already with six time-delays, the effect of a quintessence cosmology on the estimation of Ho is therefore not 
negligible. A qualitative idea of how the method can perform in the future can be gained assuming that the same six 
time delays can be estimated with only a 10% error on the variable /. In this case, we would get not only a better 
estimation of Ho but also a substantial removal of the degeneracy with respect to fl m and m. This exercise illustrates 
the potentiality of the method towards a detection of quintessence and a distinction from a pure cosmological constant. 
The results of the simulation for a = are summarized in Table IV and in Figure 6. 



Table IV 





Max 


1(7 


2a 


3ct 


H (km/s/Mpc) 


62 


59 < Ho < 66 


56<H < 71 


55 < H < 76 







n m < o.i3 


fl m < 0.37 


n m < 0.59 


to 


3 


to > 1.62 


m > 0.44 


TO > 0.10 



IV. CONCLUSIONS 



If 70% or so of the total matter content is filled by a new component with negative pressure and weak clustering, 
all the classical deep cosmological probes are affected in some way. Here we addressed the question of how the 
Dyer-Roeder distance changes when this new component, quintessence, is taken into account. We have shown that, 
particularly in the case of empty beam, the effect of the quintessence is to move the estimate of Ho to lower values 
with respect to two standard models, and to increase the spread of the likelihood distribution. This is the first time, 
to our knowledge, that a full likelihood analysis of the almost entire set of time-delays available is performed. As a 
byproduct of our analysis, we produced fits for a large range of values of £l m and to accurate to within 5%. 

The future prospects seem interesting for the time-delay method: with a not unrealistic increase in accuracy (or in 
the number of time delays), the quintessence could be detected and distinguished from a pure cosmological constant, 
thanks to the deepness to which lensing effects are observable. 

An obvious improvement of our analysis, currently underway, is to investigate the dependence on a ^ 0, 1 and on 
curvature, producing a fully marginalized likelihood for Ho- 
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FIG. 2. The solid lines are the numerical solutions for Q m = 0.3 and, from top to bottom, for m = 0.3; 1; 1.5; 2. The dashed 
lines are our fits. 
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FIG. 3. Top panel: comparison of the measured time delays / of Table II (points with error bars) with the theoretical time 
delays (lines) for an empty beam and H — 71 km/s/Mpc. The theoretical parameters are: Q m = 0.1 and m = (solid line) 
and m = 3 (dotted line); Q m = 0.9 and m = (short-dashed line) and m = 3 (long-dashed line). Bottom panel: the same as 
above but for a filled beam and Ho = 64 km/s/Mpc. 
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FIG. 5. Top panel: comparison of the fully marginalized likelihood L(Hq) (solid line) for a = with the likelihood obtained 
fixing the cosmology to two standard cases: pure cosmological constant ( fi m = 0.3 , m = 0, short-dashed line) and no 
quintessence (Q m = 1 , m = 3, long-dashed line). Notice the significant shift. Bottom panel: same as above but for a = 1. 
Here the shift is reduced. 
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